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We study the effect of strong electron-electron interactions on the persistent cur- 
rent in a multichannel ring with the aid of the Bethe ansatz equation discussed by 
Sutherland p] and Tsvelik [||. The interaction acts differently on a single and mul- 
tichannel ring. If, in the first case, the current is suppressed, then in the second 
case the current may be enhanced. With the interaction the current is proportional 
to square of the number of channels. The maximum effect corresponds to medium 
values of the interaction, where U ~ Vp. The new phenomenon of the fractional p/q 
Aharonov-Bohm effect appears for interacting electrons on the multichannel ring. 



In 1959 Aharonov and Bohm proposed an experiment, where they described the 
periodic dependence of an electron beam difraction pattern on the magnetic flux of the 
tiny solenoid penetrating the beam. One year later such single flux quantum periodic de- 
pendence was observed. There is also a single flux quantum periodic dependence for a 
persistent current in a small metallic ring, penetrated by a magnetic flux, similar to that 
in the Aharonov-Bohm experiment, which is also known as Aharonov-Bohm (AB) effect [|J, 

i- 0- 

Now, let us discuss the role of an electron-electron interaction in the AB effect. First 
of all, we consider the single channel ring with spinless fermions in a tight-binding model, 
where there is only a repulsion between electrons located on next - neighbor sites. The 
numerical simulation indicates that there is no much difference between short and long- 



range interaction ||11|| . It is intuitively clear that the interaction induces a repulsion between 
levels. The latter causes the Fermi velocity Vp and hence the persistent current to be 
reduced. However, let us note that the continuum and lattice model give different results 
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11 1 . For the continuum single channel model the two-particle interaction does not change 



the energy-flux dependence [[II]], which is in a contrast with the lattice models |[PT |. 

In weak coupling this effect may be shown with the aid of a bosonisation in the tight- 
binding model, where there is the next neighbor repulsion W. As a result, in the expression 
for the persistent current for free fermions: 

Jp~-¥jr(f-p) (CI) 

we must only change the Fermi velocity 



/ W 

V F ^U F = V F Jl-— (0.2) 
where the parameter p is an integer or half-integer number( we use the units e — h — c—1). 



The current depends on the parity of the total number of fermions iV |T7|], [18], [|T^]. However, 



in the continuum model the Fermi velocity is unchanged JT9[ . 

Thus, at small W/Vf the suppression of the current is almost linear in W. However, 
with increasing W, the suppression is saturated and the current amplitude takes a limiting 
value depending on the number of fermions on the ring. Applying the Bethe ansatz to the 



limit W/Vf — > oo, we have the following exact expression for the current [ l7| : 



4vr sin 2 vr 
J P = - T —^sm T (f-p) (0.3) 

The impurities on a single channel ring induce a repulsion between levels, therefore 
one may expect the suppression of the persistent current. However, the interaction may 
effectively renormalize the localization length, which slightly increases, i.e. there might be 
a little enhancement ( see, for comparison, Ref. [[UJ |LlJ pU| [fHp. 

Thus for a single channel ring 1) there occurs a suppression of the persistent current; 2) 
with the disorder the suppression is slightly reduced or enhanced. 3) there exists a parity 



effect: a difference in the current behavior for an even and an odd number of particles ||17|| , 

However, the situation is drastically changed when there are two or many channels or 
even if we take into account the spin quantum numbers. To show this let us consider a ring 
consisting of quantum dots. The system has effectively two channels for up-spin and down- 
spin particles, and is described by Hubbard model with t as a hopping element between 
dots, and U as a charging energy U = e 2 /C, where C is the capacitance or radius of a single 
dot. Then, the Hamiltonian reads 

H = -tH C i,«C i+ x,„ + t/£ Ni } N i{ (0.4) 

2,(7 i 

where U is a characteristic energy of a level splitting on a single dot. 

One may strightforwardly estimate the persistent current in the weak coupling, i.e. when 



U/Vf « 1 via bosonisation |pffi| . Doing so, we see that in this case the spin and charge 
degrees of freedom are decouple, i.e. there are spinons and holons and there is no suppression 
of the persistent current. The effective Fermi velocity of charge degrees of freedom (holons), 
characterising the current, does not decrease, like that for spinless fermions, but, instead, 
increases 



V f ^Ucf = V f Ji + ¥- (0.5) 
In contrast with the holon Fermi velocity, the Fermi velocity of spinons decreases Vp — > 



The persistent current is carried by holons, having a charge but no spin. Note that 
to create a single holon we have to have two electrons, as follows from bosonisation. For 
examle, a single holon is equal to an electron with up-spin f plus an electron with down- 
spin I, resulting in a charge without a spin ( a fluctuation of the charge density, only). 



Therefore, the Fermi momentum of holons is 2kp. Therefore, the Aharonov-Bohm effect 
or persistent current is characterized by half-flux quantum periodicity. Thus, the reason of 
such periodicity is a decoupling of the spin and charge degrees of freedom. 

To get the solution in the strong-coupling limit, we expand the Bethe equation in the 
parameter a = Vp/U = tN/UL << 1, which is equivalent to a large ratio U/t or a low 
density of electrons N/L << 1. In the zeroth approximation the spectrum is ||17|| : 

K n =f(n + f + ^), (0.6) 

which is different from that for spinless fermions, which is K n = ^(n + /), by an additional 
statistical flux $si = J2a Ja/N, where the J a are the quantum numbers of spinons. The 
statistical flux is changed by fractions l/N, where I are arbitrary integer numbers. This 
means that the persistent current is a l/N periodic function |2^] and described by the 
formula: 

Jp~~(f-P) (0-7) 
where p = l/N. In other words, the current displays the fractional jr Aharonov-Bohm effect 



||23|| , |[24|| , ||25|| . Note, that the holes have multivalued wave function, which may be related 
to fractional statistics. 

When N — > oo the amplitude of the fractional l/N oscillation of the current vanishes, 
i.e. the fractional effect disappears. Then, the main periodicity will come from the first 
correction in the parameter a. In spite of the fact that the spin and charge degrees of 
freedom are decoupled the spinons create a fictitious field (j)f ic , in which the holons move, 

/ic = ^§^ = £^ = ^. (0.8) 

n 

where K s is the total momentum of spinons. 



To estimate the dependence of the spinon energy on the fictitious flux we will use the 
Bethe ansatz equations(see, |26[| and references therein). First of all, let us consider the case 



when the spinon band is half-filled, that is, when the number of up-spin particlesiVf = M = 
N/2. 

In the expansion in the parameter a, the spinons are described effectively by a Heisenberg 
Hamiltonian p6 |. The dependence of the spinons' energy on the total momentum (fictitious 



flux) is the energy of the spin-wave excitations, so we get 

VpOL 

E s = —j— | sm0/ ic | . (0.9) 

The total energy of the Hubbard ring is the sum of the holon and the spinon energy: 

E = - 0) 2 + ^fa | sin0 | (0.10) 

where fictitious flux takes the values = 0, ... 

With a change of the flux /, at zero fictitious flux <p the holon energy (the first term 
in eq. flU.lUfl ) increases rapidly. However, if with the change of / we also change 0, exciting 
spin-waves, the contribution of the first term may vanish and the total energy decreases. 
Thus in the large N limit, we must put <f> = f and the ground state energy dependence is 
described by the second term in the eq. ( P-10|) . This dependence is a |-flux quantum periodic 
function, which is an envelope of -^-flux periodic parabolas. Therefore, the persistent current 
is described by the expression 

V F 

J p ~ — —acosf (0-11) 

Lj 

where < / < it. There is also a difference (a quarter flux quantum shift) in the behavior 
of the persistent current for even and for odd N (parity effect). 
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So, we arrive at the conclusion that the structure of the spin-wave excitation spectrum, 
i.e. the dependence on the spin wave momentum, determines the AB periodicity as well as 
the parity effect. Comparing the formula ( p.5| ) with the formulae( p.ll| ) one sees that the 
amplitude of the current is maximal at intermediate values of U, when U ~ Vp. That is, on 
the two channel ring there is an enhancement of the persistent current, which is maximal at 
intermediate values of the electron-electron interactions. 

In the case of arbitrary M we find that the envelope function of 1/N flux quantum 
periodic parabolic curves is an M/N periodic harmonic (for details, see Ref. \W^ )- Strictly 
speaking, there are M/N and single flux quantum periodic oscillations of the persistent 
current. The M/N oscillations survive even if the a increases, see, Fig.l, where we show the 
1/3 periodic oscillations of the ground state energy flux dependence ( low energy parabolic- 
cusp curve). The reason of such periodicity is a cuspoidal shape of the spin wave excitation 
spectrum, which is presented by short horisontal lines on this Figure. One sees that in this 
case the persistent current is determined by lowest parabolas and, therefore, its maximal 
value is proportional to ~ lOVp/L. 

However, when a descreases, the energy of spinons decreases and at some flux value this 
energy becomes the lowest one ( see, Fig.2A ). Then, in this region the 1/N oscillations are 
created which envelope is determined via the spinon excitation energy. Therefore, in this 
case instead of parabolic-cusp function (see, Fig.l) the ground state energy-flux dependence 
becomes a smooth function (see, Fig.2B). 

The most stricking result, however, is that the persistent current is a perfect 1/3 periodic 
function of the flux in units of elementary flux quantum ( see, Fig. 3). There are two 
different regions, where the current effectively increases and monotoneosly decreases. In 
the flux region where the current effectively increases there occurs the perfect 1/N periodic 
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oscillations, which are practically not resolved on the scale of the Figure. The current 
decreases linearly on the each 1/N flux region. The slope of this decrease in the single 1/N 
oscillation is equal to the slope in the region where the current monotoneously decreases 
(see, Fig. 3). The maximal current amplitude is equal to ~ 10Vp/L. Note, that this large 
amplitude of the current is determined via the small 1/N— oscillations, related to the N— 



particle bound state p3fl . Therefore on the Hubbard ring with the disorder the localization 
effects are weaker, since the localization length of this iV— particle bound state is in N— 
times larger. 

Generally speaking, the amplitude of the 1/N current oscillations is proportional to 
~ Vp/LN. The amplitude of M/N flux quantum periodic current oscillations is proportional 
to ~ Vpot/L. Thus, the M/N periodic current oscillations are dominant if a » 1/N ( see, 
for details Ref. j27|). 

It is obvious that in the more complicated double chain ring, the fractional M/N periodic 
oscillations can occur, provided that there is an effective separation of the " spin and charge" 
degrees of freedom. The " spinons" , whose number equals M, are associated with the relative 
fluctuations of the charge density and the "ho Ion" — with the total fluctuation of the charge 
density in these two channels. As a result, the total momentum of "spinons" is a fictitious 
flux for the motion of the "holons" . This fictitious flux creates the fine structure of AB effect 
on the double chain ring. 

There is a generalization of the results obtained to a multichannel ring. So, let us assume 
that there are p-channels, where there arises a relative splitting of degrees of freedom into 
holons and p — 1 types of "spinons". The problem may be studied with the aid of the 
Bethe ansatz for the SU(N) * SU(M) model. The Bethe ansatz equation for this model has 
been discussed by Tsvelik [[J. We just consider the SU(p) version of this model, which for 
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the repulsive interaction was first studied by Sutherland U and insert the flux via twisted 
boundary condition. For an attractive interaction this model has been studied by Takahashi 
p8| . The model allows to find the exact analytical solution in the limit of the strong 
coupling Vf/U « 1, which is a generalization of the discussed solution for the Hubbard 
ring presented above. Generally speaking, in the expansion with the parameter a the Bethe 
ansatz equations of the SU(p) model are decoupled into a system of p — 1 dimensionless 
equations, which describe p — 1 fermion reservoirs of " spinons" and a single equation which 
describes the charge degrees of freedom. 

There occurs (p — 1) fictitious fluxes fa associated with p — 1 types of "spinons", which 
are relative fluctuations of the charge density between channels. In other words, one may 
say that there are Mi, ..,M p _i spinons, located on the p — 1 channels. The "holons", the 
number of which is equal to the total number of particles N move in the fictitious field 
created by the p — 1 sorts of spinons. 

The total energy of the multichannel ring in magnetic field extracted from the Bethe 
ansatz solution has the form: 

E = T (/ " 01 " 02 " - " 0p ~ l)2 " T aF(01 ' 02 ' (ai2) 

where the fictitious fluxes take the values: fa = 2 ^ i , and Ki = 0, ±1, ±2, .... The function 
.F(0i, 02) •••> fa-i) describes the energy of p — 1 reservoirs of spinless fermions(" spinons"), 
interacting with each other with a dimesionless coupling. The form of F{fa, fa, <f> p -i) 
depends strongly on the distribution of particles between channels. Therefore, there are p—1 
different Fermi momenta associated with the numbers Mi, .., M p _i where, for example, Mi is 
the number of particles in the i— th channel. For this reason, the multidimensional function 
F(fa,4>2, may have cuspoidal minima when the flux value fa commensurates with 
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the appropriate Fermi momentum, i.e. with kf . 

The presented expression (p.!2|) allows us to get the following results for the multichannel 
ring: analogous to the case of the two channel ring, there is an enhancement of the persistent 
current for an intermediate value of the electron-electron interaction; 2) there are (p — 1) 
fractional periodicities Mi/N, related to the particle distribution between channels. 

The other important point is the dependence of the persistent current on the number 
of channels. Note that Vp = tN/L is related to the total number of particles N on the 
ring and therefore the Fermi velocity in a single channel vf is equal to Vp/p. That is, the 
amplitude of the persistent curent is proportional to v 2 F p 2 / (LU) and increases as p 2 with the 
number of channels. Therefore, with an increase of the number of channels one may expect 
a very strong enhancement of the persistent current. The disorder, however, suppresses 
this enhancement, but the amplitude of the persistent current may still remain very large 
and may be compared with the current of the free electrons. Whether or not this effect 



can explain the discrepancy between experiments [|13j and [JT2J] and existing theories is still 
unclear. 

However, one may suggest that the large amplitude of the current on a single metallic ring 
I2| is due to an effect described in the present paper ( due to an large number of channels). 



The small value of the current in the Levy et al experiment (~ 10~ 2 vf/L) is due to an average 
over an ensemble of rings with different number of particles. Such averaging explains well 
both the half- flux quantum periodicity as well as the suppression of the persistent current. 
In the latter experiment the current on the different rings has different orientations resulting 
in a compensating effect. Since the distribution of particles over the channels is, probably, 
well controllable, it should be a realizable experiment to see the predicted fractional periods. 
Recently we became aware of the work of Avishai and Berkovits ||29||, where they also 
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have seen in numerical simulations of small disordered rings with 2 and 3 channels an 
enhancement of the persistent current. 
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Figure Captions 

Fig.l The behavior of the ground state energy as a function of flux / ( the lower 
parabolic-cusp curve) and the spinon energy ( indicated by the short horisontal lines) for 
303 electrons with the 101 up spin electrons at the values L = 1000 and U — 10 in the 
region within the half of fundamental flux quantum. The energy is expressed in the units 
tlO 2 . The zero energy corresponds to — 520. Ot. 

Fig. 2 A) The behavior of the holon energy ( the parabolic-cusp curve) and the spinon 
energy, indicated by short horisontal lines as a function of flux / ($) for 303 electrons 
(M = 101) at the values L = 20000 and U = 10 in the region within the half of fundamental 
flux quantum. The zero energy corresponds to — 605.771443t. The energy is expressed in 
the units tlO 7 . B) The behavior of the ground state energy as a function of flux /. 

Fig. 3 The behavior of the persistent current as a function of flux / for 303 electrons 
at the values L = 20000 and U — 10 in the region of flux within the half fundamental 
flux quantum. The 101 particles have up-spin. The maximal current amplitude is equal to 
~ 10Vf/L. On the slope where the current increases, there are perfect 1/303 flux quantum 
periodic oscillations. 
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